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ABSTRACT 

Future microwave sky surveys will have the sensitivity to detect the kinematic Sunyaev-Zeldovich 
signal from moving galaxy clusters, thus providing a direct measurement of their line-of-sight peculiar 
velocity. We show that cluster peculiar velocity statistics applied to foreseeable surveys will put sig- 
nificant constraints on fundamental cosmological parameters. We consider three statistical quantities 
that can be constructed from a cluster peculiar velocity catalog: the probability density function, the 
mean pairwise streaming velocity, and the pairwise velocity dispersion. These quantities are applied 
to an envisioned data set which measures line-of-sight cluster velocities with normal errors of 100 km/s 
for all clusters with masses larger than 10 14 solar masses over a sky area of up to 5000 square degrees. 
A simple Fisher matrix analysis of this survey shows that the normalization of the matter power 
spectrum and the dark energy equation of state can be constrained to better than 10 percent, and 
the Hubble constant and the primordial power spectrum index can be constrained to a few percent, 
independent of any other cosmological observations. We also find that the current constraint on the 
power spectrum normalization can be improved by more than a factor of two using data from a 400 
square degree survey and WMAP third-year priors. We also show how the constraints on cosmological 
parameters changes if cluster velocities are measured with normal errors of 300 km/s. 
Subject headings: Cosmology: theory, Cosmology: cosmological parameters, galaxies: clusters, statis- 
tics 



1. INTRODUCTION 

One of the great triumphs of cosmology has been 
a precise determination of the parameters defin- 
ing the standard cosmological model, primarily via 
the cosmic microwave background fluctuations com- 
bined with the large-scale distribution of galaxies 
and the distance-redshift relation of distant super- 
novae. Future high-precision probes of cosmological 
parameters using different sources of data provide the 
possibility of tighter constraints on the cosmological 
parameters, along with strong consistency checks of 
the cosmological model. Here we consider one such 
potential probe, the peculiar velocities of galaxy clus- 
ters. U pcoming microwave sk y s urveys like Planck , 
ACT (Kosowsky et al. 2006t iFowler et al. 2005t 



mass limit. Massive clusters are the largest gravitation- 
ally bound objects in the universe, and their distribu- 
tion and velocities probe the mildly nonlinear regime of 
structure formation. For given initial perturbations in 
the early universe, peculiar velocities at recent epochs 
depend on the expansion history of the universe. Result- 
ing constraints on cosmological parameters will comple- 
ment those from microw ave background measurements 
(e.g. ( Spergel et al 2006h ^ and from supernova searches 
(e.g. ( Astier et al. 2006f h. Cluster number counts from 



iKosowskv 20031 ) and SPT (jRuhl et al 2004h will detect 
thousands of galaxy clusters via their thermal Sunyaev- 
Zeldovich distortion s of the microwave back ground 
frequency spectrum (|Sunyaev fc Zeldovich 1980l ). With 
sufficient sensitivity, angular resolution, and control 
of systematic errors, microwave surveys also have 
the potential to measure directly the line-of-sight 
peculiar velocity of clusters via their kinematic Sunyaev- 
Zeldovich signature, a near-blackbody distortion of 
the microwave background proportional to the optical 
depth and velocity of the cluster. Since the SZ effect is 
effectively independent of cluster redshift, such a survey 
can probe the cosmic velocity field over a cosmological 
volume of the universe. 

In practice, future data sets will give sky positions, 
redshifts, and the line-of-sight peculiar velocity compo- 
nent for massive galaxy clusters above a certain SZ flux 
threshold, which will roughly correlate with a cluster 
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their thermal SZ signature also show this complementar- 
ity, but they are prone to bias from systematic errors in 
cluster mass determination (Francis, Bean & Kosowsky 
2005). Cluster peculiar velocities are much less sensi- 
tive to this systematic error, and so may provide a more 
promising route to cosmological constraints, even though 
the kinematic SZ signal is much smaller than the thermal 
SZ signal. 

Other recent studies in this area (|Peel 20061: 
iMonteagudo et al 20061 : iPeel fc Knox 2002ft laid out the 
theoretical framework for studying the correlation func- 
tion of the cosmic velocity field through the cluster kSZ 
effect while DeDeo, Spergel & Trac 2005 studied cross- 
correlations of kSZ signal with galaxy counts. This paper 
explores how well various velocity statistics for a galaxy 
cluster velocity survey can constrain cosmology. 

2. THEORY 

We will consider three different velocity statistics: the 
probability density function of the line-of-sight compo- 
nent of peculiar velocities f(v), and the mean pairwise 
streaming velocity My (r) and pairwise velocity disper- 
sion <Jij(r) for two clusters separated by a distance r. 
These are formed by averaging over all clusters or clus- 
ter pairs in a given peculiar velocity survey. In the 
halo model, these quantities can be written as the sum 
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of the contribution from one-halo and two- halo terms. 
Since we are interested in very massive clusters, the 
one-halo term can be neglected. When compared with 
the VIRGO numerical simulation of struc ture formation 
(jEvrard et al. 2002t IColberg et al. 2000( 1. the two-halo 
term approximates the probability density function fairly 
well, and for Vij(r) and Uijir) it agrees well for separa- 
tions r larger than 20 Mpc and 30 Mpc, respectively. 

Analytic approximations for all three of these statistics 
in dark matter-dominated cosmological models give their 
dependence on parameters describing the background 
cosmology. The probability density function of the line- 
of-sight pec uliar velocity componen t at some redshift z 
is given by (|Sheth fc Diaferio 20011) 



/(«,*) = 



J dm mn(m\6)p(v\m, 5) 
J dm mn(m\S) 



(1) 



where n(m\8)dm is the number density of halos that have 
mass between m and m + dm in a region with overdensity 
6, and p(v\m, 5) is the probability that halos of mass m 
in a local overdensity 5 move with line-of-sight velocity 
v; the dependence of these quantities on redshift is left 
implicit. Assuming Gaussian initial conditions, we can 
write n(m\5) in terms of the total number density n(m) 
and a large-scale bias factor b(m), 

n(m\6) » [1 + b(m)8] n(m). (2) 

Define moments of the i nitial mass distribu tion with 
power spectrum P(k) by (jBardeen et al. 19861 ) 
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when smoothed on the scale R(m) = (3m/47rpo) with 
the top-hat filter W(x) — 3[sin(x) — x cos(x)]/x 3 , and po 
the present mean matter density. Th en the bias factor 
can be written as (jSheth et al 2001a( l 

a 2 (m) 



b(m) = 1 + 



S 2 ■ 

u crit 



a 2 {m)S clit D a 



(4) 



where D a is the linear growth factor at scale factor a, 
normalized to 1 today, and the critical overdensity <5 cr it ~ 
1.686. For the number density of halos n(m ), we use the 
Jenkins mass function (|Jenkins et al. 20 01) 



dn 
dm 
x exp 



(m,z) = 0.315 



Po dhi<Jo(m) 
m 2 d\nm 



|0.61 -ln((T (m)D„)| 



(5) 



This mass function is a fit to numerical simulations of 
cold dark matter gravitational clustering. 

Along with n(m\5) just define d, we also need p(v\m, 5) , 
which we assume to be gaussian (|Sheth fc Diaferio 20011) . 



p(v/m, 5) 



exp 



V2^(m)/3 



(6) 



with the three-dimensional velocity dispersion smoothed 
over a length scale R( m) given by 
(jHamana et al. 2003( 1 



<r v (m, a) = [1 + 5(R local )} 2f * {Rl - l) aH(a)D a 



d\nD a 
din a 



x 1- 



cr^(m) 



cr 2 (m)a_ 1 (m) 



1/2 



o--i(m) (7) 



Here H(a) is the Hubble parameter as a function 
of scale factor, i?i oca i is a smoothing scale with 
which the local background density 5 is defined, and 
n(R ]nra] ) =0.6a 2 (R ]nr ^ )/a 2 (10 Mpc/h) 
(She th fc Dia ferio 2003). Following 
(jHamana et al. 2003( 1. -Riocai is obtained empiri- 
cally using N-body simulations via the condition 
co(-Riocai) = 0.5(1 + z)~ - 5 . We have now completely 
specified f(v, z) in terms of quantities related to the 
background cosmology. 

The mean relative peculiar velocity v%2 (r) for all pairs 
of halos at comoving separation r can be related to the 
linear two-point correlation function £} m (r, a) for dark 
matter us ing large-scale bias and the pair conservation 
equation (|Sheth et al 2001a[) : 

v 12 (r) = --aH(a) ^ 1 + ^Un (r>o) W 

where 6haio = J dm mn(m)b(m) / J dmn(m) is the aver- 
age halo bias factor and f m (r, a) is the two-point corre- 
lation function averaged over a ball of radius r. Under 
linear evolution, the correlation function simply scales as 



e n (r,a)=D 2 a e n (r,a ) 



(9) 



Similarly, a model for the pairwise velocity dispersion 
(7i2(r), in the large-scal e limit and neglec ting velocity 
correlations, is given by dSheth et al. 20011) 



4« = 2 



(a 2 ,(R,a))+e n (r,a)(b<J 2 (R,a)) 



(10) 



where (• ■ •) represents an ensemble average over all ha- 
los lying within a given mass range. Note the above 
expressions for v\2(r) and g f 7 -(r) are for three d imen- 
sional velocity distributions. (jFerreira et al. 1999) give 
one technique for constructing an estimator of v±2 (r) us- 
ing only linc-of-sight velocities; estimators for both Vi2(r) 
and a 2 ^ (r) will be presented elsewhere (Bhattacharya et 
al. in prep.). 

These analytic approximations to the cluster veloc- 
ity probability distribution function, the cluster pairwise 
mean relative velocity, and the pairwise mean velocity 
dispersion explicitly display their cosmological depen- 
dences through the Hubble parameter H(a), the growth 
factor D a , the mean mass density po, and the power spec- 
trum P(k). Note that they depend on the background 
cosmology in significantly different ways. 

3. ESTIMATING PARAMETER CONSTRAINTS 

To understand the dependence of the distribution of 
velocities on the underlying cosmology, we perform a sim- 
ple Fisher matrix analysis for the three velocity statistics 
described in the previous Section. Assuming a flat Uni- 
verse, the set of cosmological parameters p on which the 
velocity field depends are the normalization of the matter 
power spectrum, conventionally parameterized by as, the 
power law index of the primordial power spectrum ns, 
the total matter density il m , the dark energy equation 
of state w and the Hubble parameter h. We consider 
a fiducial model with og = 0.9, ns = 1, Q m = 0.3, 
w = —1, h = 0.7. Upcoming data sets will provide clus- 
ter number counts and line-of-sight velocities as a func- 
tion of redshift; we assume here that all clusters with 



Cosmological Constraints from Cluster Velocities 



3 



Parameter 


5000 deg 2 


1000 deg 2 


400 deg 2 


Parameter 


/(«) 


iHj(r) 


cry(r) 


a 8 [5.6%] 
n s [1.5%] 
Q m [6.7%] 
w [10%] 
h [4.3%] 


7.2% {2.1%} 
2.2% {1.2%} 
7.9% {4.7%} 
9.2% {6.0%} 
18.9% {4.0%} 


16.0% {2.4%} 
4.9% {1.4%} 
17.9% {6.0%} 
20.5% {8.6%} 
42.4% {4.1%} 


25.6% {2.4%} 
7.7% {1.4%} 
28.2% {6.3%} 
31.6% {9.3%} 
67.0% {4.0%} 


(78 [5.6%] 
n s [1.5%] 

n m [6.7%] 

w [10%] 
h [4.3%] 


7.3% {2.2%} 

2.9%{1.4%} 

10.5%{6.0%} 

12.3%{7.0%} 

23.0%{4.0%} 


61%{5.0%} 

3.3%{1.1%} 

25%{5.5%} 

24.7%{8.3%} 

56.2%{4.0%} 


36.8%{3.7%} 
28.1%{1.4%} 
21.1%{6.1%} 
11.1%{7.8%} 
5.6%{2.9%} 



TABLE 1 

ler errors in cosmological parameters from f(v) for clusters with 
masses M > 10 14 Mq, measured line-of-sight velocity error 
<r„ = 100 km/s, and z < 1.5, for three survey areas. Percent 
errors are given in terms of the fiducial values. Errors in braces 
include a Gaussian prior from current measurements, given in 
square brackets. 



Parameter 


5000 deg 2 


1000 deg 2 


400 deg 2 


o- 8 [5.6%] 
n s [1.5%] 
O m [6.7%] 
w [10%] 
h [4.3%] 


34.8% {3.5%} 
1.87% {0.8%} 
18.2% {4.6%} 
16.7% {7.6%} 
42.4% {3.6%} 


174.0% {3.9%} 
9.4% {1.3%} 
92.0% {6.3%} 
84.3% {9.64%} 
212.0% {3.75%} 


435.0% {4.7%} 
23.5% {1.4%} 
230.7% {6.5%} 
211.0% {9.8%} 
531.0% {4.0%} 




TABLE 2 
Same as in Table 1, but for Vij(r). 


Parameter 


5000 deg 2 


1000 deg 2 


400 deg 2 


o-8 [5.6%] 
n s [1.5%] 
n m [6.7%] 
w [10%] 
h [4.3%] 


20.2% {3.0%} 
14.8% {1.4%} 
14.1% {5.6%} 
8.3% {5.6%} 
3.2% {2.3%} 


102.0% {3.5%} 
74.7% {1.4%} 
71.2% {6.2%} 
41.6% {9.4%} 
16.2% {3.8%} 


255.0% {3.7%} 
187.0% {1.4%} 
176.0% {6.3%} 
104.0% {9.8%} 
40.4% {3.9%} 



TABLE 3 
Same as in Table 1, but for aij(r). 



M > 10 14 M Q are detected, that redshifts are measured 
with zero error for z < 1.5 and that peculiar velocity 
is measured with a standard error of 100 km/sec, the 
limit from internal cluster motions (Nagai, Kravtsov & 
Kosowsky 2003). This is likely a re alistic mass limit fo r 
high-sensitivity microwave surveys (|Sehgal et al. 20"07h ■ 
The velocity error represents the fundamental limit from 
intrinsic cluster properties; measurements to this pre- 
cision will likely require measurements of cluster tem- 
perature (Sehgal, Kosowsky & Holder 2005) and de- 
tailed understanding of the point sources and other sig- 
nals which can co n taminate the kinematic S Z signal 
dCoble et al. 20061: iHoIder fc Kosowsky 20041) . How- 
ever (jDiaferio et al. 2 005) have extracted velocities 
from cluster simulations using multiple microwave and 
X-ray bands finding a robust lower limit of 200 km/s. 
We also show how constraints on different parameters 
degrades for a larger standard error of 300 km/s. Spec- 
troscopic redshifts will be challenging to obtain for clus- 
ters beyond z = 1.5, so we consider clusters only to this 
distance. We approximate the cosmic variance error for 
each observables to be proportional to (r/L) 3 , where r is 
the separation between two halos and L is the size of a 
particular survey. A detail study of cosmic variance will 
be discussed elsewhere. 

The measured observables are the value of f(v) in red- 
shift and velocity bins, and the values of (r) and aij (r) 
in redshift and distance bins; the variance in each bin 



TABLE 4 

la errors in cosmological parameters obtained from the three 
velocity statistics for 5000 deg 2 with measured line-of-sight 
velocity error a v = 300 km/s. 

depends on Poisson errors, cosmic variance, and veloc- 
ity errors. The redshift bins Az = 0.2 are chosen to be 
narrow enough to capture statistically significant redshift 
evolution. Then assuming a normal error distribution in 
each bin, the Fisher information matrix for each statistic 
is simply (Jungma n et al. 1996[ ) 



^ d(j){xk,zi) 1 d<t>(xk,zi) 

*ap = 2^ 



8p a 



(11) 



where the values zi define the redshift bins, the values 
Xk define the bins in either velocity (for /) or distance 
(for Vij and <Jij), <j> is either f(v, z), Vij(r, z) or <xy (r, z), 
is the variance of the function <j> in each bin, and the 
partial derivatives are evaluated for the fiducial values of 
the cosmological parameters. The inverse of the Fisher 
matrix then gives an estimate for the variances of each 
cosmological parameter, marginalized over the values of 
the other parameters, along its diagonal, and correlations 
between the parameters in its non-diagonal elements. 

We have looked at parameter constraints for cluster 
velocity surveys of three different sizes: sky areas of 400, 
1000, and 5000 square degrees. We calculate the pa- 
rameter errors for each of the three velocity statistics 
individually, and also when imposing prior constraints 
on each parameter taken from the WMAP 3-year data 
(jSpergel et al 20061 ) plus curr ent large-scale str ucture 
and supernova measurements (jAstier et al. 2 006). Re- 
sults are summarized in Tables 1,2, and 3 for each of the 
three velocity statistics. Table 4 shows the constraints 
for a 5000 deg 2 survey when the velocity error is in- 
creased to 300 km/s. Comparing the results for 300 
km/s with those for 100 km/s, we find that for f(v) the 
errors on parameters changes by 3-4 % while for Uy(r) 
and crij(r) this degradation is more pronounced. How- 
ever note that most of the constraints on cosmological 
parameters comes from f(v), so the overall errors on cos- 
mological parameters degrades by 3-4%. 

It is clear from these results that a number of cos- 
mological parameters can be well constrained by large 
cluster peculiar velocity surveys. In particular, both 
f(y) and Vij can provide strong constraints on ns] for 
a survey covering 5000 square degrees, the constraints of 
2.2% and 1.9% are close to the current precision from 
microwave background measurements. Note that these 
constraints on ns are at a smaller physical scale than 
the corresponding constraints from primary microwave 
anisotropies, and a mismatch in these two high-precision 
measurements indicates a departure from a pure power 
law of the primordial density perturbations. The Hub- 
ble parameter can be measured using tr^ to 3.2%, better 
than current constraints. All three statistics give con- 
straints on w to around 10%, competitive with current 
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measurements. Finally, f(v) gives a measurement of a a 
to 7.2%, competitive with current limits. These limits 
are for each parameter, marginalized over the others. 

4. CONCLUSION AND PROSPECTS 

Galaxy cluster peculiar velocities, measured directly 
via the kinematic Sunyaev-Zeldovich effect, represent a 
new route to precision cosmological constraints. We have 
shown here that if cluster line-of-sight velocities can be 
measured with errors of 100 km/sec, the resulting con- 
straints on several cosmological parameters will be com- 
parable with all current techniques. Such measurements 
could be important for their constraints on particular 
parameters, but are likely more valuable as consistency 
checks on the standard cosmological model. Multiple 
measurements of each cosmological parameter using in- 
dependent methods provides our only way to determine 
whether our universe is actually described by the simple 
models spanned by the standard cosmological parameter 
space. The awkward appearance of dark energy on the 
cosmological stage makes these cross-checks all the more 
imperative. 

A potentially even more important advantage to pe- 
culiar velocities is the control of systematic errors. It 
has been appreciated for some time that the comoving 
number density of clusters above a certain mass limit 
as a function of redshift depends sensitively on the un- 
derlying cosmology: a small increase in the cosmological 
growth factor leads to a large increase in the number of 
clusters (Viana & Liddle 1996; Eke et al 1996; Barbosa 
et al 1996; Bahcall, Fan & Cen 1997). Blind Sunyaev- 
Zeldovich surveys will detect all clusters above a cer- 
tain SZ-distortion threshold in a given direction of the 
sky, so their signal also depends sensitively on cosmolog- 
ical parameters. But the cluster SZ selection function 
is not equivalent to a simple mass cutoff, and a small 
systematic error in understanding this selection function 
can result in substantial systematic errors in cosmologi- 
cal constraints (Francis, Bean & Kosowsky 2005; Holder 
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et al 2001). In contrast, the velocity statistics considered 
here have little dependence on the cluster selection func- 
tion, as we will show explicitly elsewhere; 20% variations 
in the mass selection function give only a few percent 
change in the velocity statistics. 

The observational prerequisite for measuring galaxy 
cluster velocities are microwave observations with multi- 
ple frequencies, arcminute angular resolution, and sen- 
sitivity of a few [iK per arcminute pixel. The Ata- 
cama Cosmology T elescope (ACT) will no minally meet 
these requirements (jKosowskv et al. 20061) ; it is sched- 
uled for deployment in Chile in the first part of 2007 
(although it will not commence three-frequency obser- 
vations until 2008). The South Pole Telescope (SPT) 
also could make sufficiently sensitive maps, although its 
nominal initial surv ey will be substa ntially wider and 
shallower than ACT (|Ruhl et al 20041) . Within the next 
few years, we will have the raw experimental capabil- 
ity to construct cluster peculiar velocity catalogs. The 
extent to which this program will yield the kinds of cos- 
mological constraints discussed here will then depend on 
our ability to extract the small velocity signal from the 
mix of other larger signals at multiple wavebands, in- 
cluding the clusters' thermal SZ distortion, infrared and 
radio point sources (Coppin et al 2006; Knox, Holder & 
Church 2004), and gravita tional lensing by the cluster 
(|Holder fc Kosows ky 2004] ). The remarkable strides in 
observational cosmology over the past ten years suggest 
that optimism is justified. 
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